A Testbench for the Nested Dipole Hypothesis of Kosterlitz 
and Thouless 



A.AlastueyQ 

Ecole Normale Superieure de Lyon,Laboratoire de Physique, Unite de Recherche Associee 
1325 au Centre National de la Recherche Scientifique, 46,allee d'ltalie, 69364 Lyon Cedex 
07, France 

P.J.Forresteif| 

Department of Mathematics, University of Melbourne, Par kville, Victoria 3052, Australia 



We consider the two-dimensional one-component plasma without a background and 
confined to a half plane near a metal wall. The particles are also subjected to 
an external potential acting perpendicular to the wall with an inverse power law 
Boltzmann factor. The model has a known solvable isotherm which exhibits a 
Kosterlitz-Thouless type transition from a conductive to an insulator phase as the 
power law is varied. This allows predictions of theoretical methods of analysing 
the Kosterlitz-Thouless transition to be compared with the exact solution. In par- 
ticular we calculate the asymptotic density profile by resumming its low fugacity 
expansion near the zero-density critical coupling in the insulator phase, and solving 
a mean-field equation deduced from the first BGY equation. Agreement with the 
exact solution is obtained. As the former calculation makes essential use of the 
nested dipole hypothesis of Kosterlitz and Thouless, the validity of this hypothesis 
is explicitly verified. 
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1. INTRODUCTION 



The two-dimensional Coulomb gas refers to a neutral system of charged particles confined 
to a plane. The two species have opposite charge (magnitude q say) and interact via the 
laws of two-dimensional electrostatics (logarithmic potential) To stop collapse between 
oppositely charged particles at low temperature, due to the singular behaviour of the 
logarithmic potential at the origin, a hard core or similar short range regularisation is 
also required. For low densities and high temperature the two-dimensional Coulomb gas 
forms a conductive phase in which the positive and negative charges are dissociated and 
can screen a long- wavelength external charge density. In contrast, for low densities and low 
temperature, the system forms a dipole phase in which the positive and negative particles 
pair together . Perfect screening of a long-wavelength external charge will no longer occur. 

Intricate structures of nested dipoles were hypothesized by Kosterlitz and Thouless 
[1] for the dominant configuration contributing to the polarization as the transition point 
is approached from the dipole phase. On the basis of this remarkable physical insight , 
an iterated mean- field theory was formulated and quantities of physical interest thereby 
calculated in the vicinity of the critical point. This so called Kosterlitz-Thouless transition 
between the conductive and dielectric states occurs at the coupling T = 4 (r := q 2 /kT)'m 
the zero-density limit. 

The nested dipole hypothesis and the iterated mean-field equations of Kosterlitz and 
Thouless were recently put on a firmer footing by Alastuey and Cornu [2], who made a low 
fugacity (£) analysis of the charge-charge correlation function and the dielectric constant 
e for T — > 4 + . At order ( 4 it was proved that the configurations giving the leading 
order contribution to 1/e are the nested dipoles hypothesized by Kosterlitz and Thouless. 
Assuming this to be true at all orders, the low fugacity series could be resummed, and 
the iterated mean field equations of Kosterlitz and Thouless derived exactly. 

The pairing transition from a conductive to a dielectric phase is not unique to the two- 
dimensional two-component Coulomb gas. One-component log-potential Coulomb gases 
also exhibit this transition, provided the neutralizing background consists of a lattice 
of oppositely charged particles [3], or there is no background and the system is in the 
vicinity of a conductive medium [4]. For the latter class of system a solvable model has 
been formulated which exhibits a pairing transition as a microscopic parameter is varied 
[5]. The model is the two-dimensional one-component plasma consisting of particles of 
positive charge only confined to a half plane in the vicinity of a metal wall, and subjected 
to a one-body external potential such that 

e -PV{y) =y -a 

which acts in the direction perpendicular to the metal wall only. It is solvable at the 
special coupling T = 2, and exhibits a pairing transition as a is varied through one. It 
is our objective herein to use the exact solution as a testbench for the predictions of 
theoretical methods of analysing the Kosterlitz-Thouless transition for this model. 

We begin in Section 2 by considering the screening properties of the system with 
respect to an infinitesimal external dipole. In Section 3 the low fugacity resummation 
technique of [2] is applied to study the density profile in the dielectric phase near criticality, 
and the density profile is further analysed using the first BGY equation. In Section 4 
comparison of the theoretical predictions with the exact results are made. Concluding 
remarks are made in Section 5. 
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2. CHARACTERIZING THE PHASE 



2.1 Definition of the model 

Consider a system of two-dimensional charges of strength q confined to a half plane y > d 
and suppose a perfect conductor occupies the half plane y < 0. For each charge of 
strength q at position (x, y) say in the system , the effect of the perfect conductor is 
to create an image charge of strength — q at position (x, — y).The electrostatic potential 
0(r, r') experienced by a test particle of charge q at r' — (x', y') due to a particle of charge 
q at r = (x, y) is then 

0(r, r') = q 2 (v c (\f— r'\) — v c (\f — r')|) (2.1a) 

where 

v c (r-P) = -log[{(x- x') 2 + (y - y') 2 } 1/2 / L] and r = (x,-y) (2.16) 

For convenience the arbitrary length scale L will be henceforth set equal to unity. As 
well as interacting via the pair potential (2.1), the particles also experience a one-body 
potential with Boltzmann factor (1.1). Since the electrostatic potential 4>(y) due to a 
background charge density qpb(y ), d <y < oo is given by 

<j>{y) = -*<l dy'(\y-y'\-(y + y'))p b (y') (2.2) 

J d 

it is straightforward to check that the one-body potential given in (1.1) can be interpreted 
as being due to a background charge density 

M) - ^ (2.3) 

However it is more convenient for our purposes below to interpret (1.1) as non-electrical 
in origin, and making no contribution to the total charge density. 

2.2 Response to an external dipole 

The conductor and dipole phases of the two-dimensional Coulomb gas can be distinguished 
by different screening properties of an infinitesimal external charge: the external charge is 
perfectly screened in the conductor phase while only a fraction 1 — 1/e is screened in the 
dipole phase. For Coulomb systems near a metal wall an external charge is automatically 
screened by its own image. We consider instead the screening of an infinitesimal dipole. 
The image of a dipole pointing perpendicular to a metal wall has the same magnitude 
and direction as the original dipole. A conductor phase should perfectly screen such an 
external dipole. 

Let us use linear response theory to give a mathematical characterization of the screen- 
ing of an external dipole. An external dipole at r := (0,yo) which is of strength p and 
perpendicular to the metal wall adds to the Hamiltonian a term 



H ex t = Po 



[ dr' ^-v c {\? - r Q \)Q{?) + [ d? ^- Vc {\? - f Q \)Q{?) 
Jv oy Jv oy 



(2.4) 



3 



where r% := (0, —yo),y~o '■= —yo and Q(r') denotes the microscopic charge density at point 
r'. The domain V is the half plane y > d. According to linear response theory the change 
in charge density at a point r, 5q(r) say, due to the external dipole is given by 

5q(r) = -fi i(H cxt Q) - (H cxt ) (Q)} (2.5) 

From (2.4) the r.h.s. of (2.5) can be written in terms of the charge-charge correlation 

S(f,r') := (Q(?)Q(r>)} - (Q(7))(Q(r>)} (2.6) 

as 

8q(r) = -Ppo I dPS(fy'){^-v c (\r'-ro\) + ^v c (\r'-r%\)\ (2.7) 
Jt> [oyo dy J 

On the other hand, our characterization of the conductor phase as perfectly screening the 
dipole says 

/ dr y5q(r) = -p (2.8) 
Jv 

Substituting (2.7) in (2.8) gives the sum rule 

(3 [ dry ( d?S(r,?){^v c (\P-f \) + ^v c (\P-f \)\ = l (2.9) 
Jv Jv [dy dy J 

to be obeyed by the system in the conductor phase. 
The sum rule (2.9) can be further simplified. Now 

d JIS r°° , d 



where 



\ dry] dr' S(r,r')—v c (\r' -f \) = dyy—F(y ,y) (2.10a) 
Jv Jv dyo Jd dyo 

/oo /*oo rco 

dx dx' dy'S(y',y;x' - x)v c (x' - x ;y' - y Q ) (2.106) 
-oo J —oo Jd 



From the convolution formula for Fourier transforms we have 

F(y ,y)= dy'S(y',y;0)v c (0;y'-y ) (2.11) 
Jd 

In performing this step we are assuming that with y fixed S(x' — x; y', y) decays sufficiently 
as a function of y' and x — x' for the integral in (2.10b) to be absolutely convergent and 
thus the order of integration to be unimportant. Next, it is a straightforward exercise to 
deduce from Poisson's equation 

V 2 v c (\f\) = -2tt 5(f) (2.12) 

that 

v c (0;y' - y ) = -n\y' - y \ (2.13) 

Substituting (2.13) in (2.11) and then substituting the resulting expression in (2.10) and 
performing the differentiation gives 

/ dry / dr'S(f,r')—v c (\r'-fh\)=n / dyy / dy' S(y', y; 0)sgn( 2/ / - y ) (2.14) 
Jv Jv dyo Jd Jd 



Performing an analogous simplification on the second term on the l.h.s. of (2.9) we thus 
deduce that (2.9) is equivalent, subject to the clustering assumption for S(x' — x;y',y) 
noted below (2.11), to the simpler sum rule 

roo roo roo 

27T/3/ dyy dy' dx' S(y' : y; x' - x) = 1 (2.15) 

Jd Jyo J—oo 

This sum rule is to be satisfied in the conductive phase of any two- dimensional Coulomb 
system separated a distance d from a metal wall (when y = this result has been obtained 
previously by Jancovici [6]). 

A remarkable property of (2.15) is that it must hold for all positions (0, yo) of the 
external dipole. Differentiating with respect to y , and changing variables x' \— > x' + x — x , 
we thus have 

/ dryS{f ,7) = 0. (2.16) 
Jv 

Hence the perfect screening of an external dipole implies the dipole moment of the internal 
screening cloud must vanish. For a phase which does not perfectly screen an external 
dipole, the sum rule (2.15) cannot hold. Rather, we would expect the l.h.s. to depend on 
y and thus the dipole moment of the internal screening cloud will be non-zero. 

2.3 Phase transition and potential drop 

In electrochemistry a fundamental quantity is the potential drop across the interface: 

roo 

A0 = 27rg/ dyyp(y) (2.17) 

J d 

where qp(y) denotes the total charge density at distance y from the interface. As previ- 
ously noted [5], the formula (2.17) exhibits a further interpretation of A0: it is directly 
proportional to the mean distance between a particle and the metal wall, or equivalently 
the mean size of the particle-image pairs. Therefore, A0 is expected to be finite in the 
insulator phase, while it should diverge in the conductor phase where the charges of the 
plasma are not paired by their own images. As seen from the integral expression (2.17), 
the finiteness of A0 is closely related to the large-distance behaviour of p(y). In the 
insulator phase, p{y) should decay as l/y 2+e (e > 0) when y — > oo, while in the con- 
ductor phase p(y) should decay typically as 1/y 2 or slower. The relation between A0 
and the asymptotics of p(y) will be studied in Section 3, by resumming the low fugacity 
expansions. 

The potential drop A<j>, or the internal dipole moment 

D(y ) :=2tt / dfyS(r ,r) (2.18) 
Jv 

may be taken as equivalent idicators for characterizing the phase of the present model. 
A relation between both quantities can be obtained by starting with the compressibility 
sum rule 

= \T dx> fV<SW;s-:iO (2.19) 

CC, Q J~oo Jd 

Taking the first moment of both sides gives 

C Hdyy^- = -T dyyT dx' f° dy' S(y,y'; x - x') (2.20) 

Jd CC, (J Jd J—oo Jd 
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Next we want to interchange the order of the y and y' integrations on the r.h.s. and 
interchange the order of integration and differentiation on the l.h.s. . From (2.17) and 
the explicit form of S(y, y';x — x') (see Section 4) we see that a necessary condition for the 
validity of both these operations is that A<j) be finite. Assuming this condition and the 
validity of the operations for the dipole phase, we obtain from (2.20), after using (2.17) 
and (2.18), the desired relationship: 



This equation can only be valid in the insulator phase, since in the conductor phase, 
the quantity on the r.h.s. of (2.20) is given by the universal value (2.15), so we have 
instead 



We stress that, in the conductor phase, the integrations over y and y' in the r.h.s. of (2.20) 
cannot be inverted. Otherwise, the corresponding integral, which also appears in the l.h.s. 
of (2.15), would vanish as a consequence of D{y) = 0. This non-absolute convergence is 
related to a slow decay of S for some configurations. At the same time, the differentiation 
with respect to ( and the integration over y in the l.h.s. of (2.20) cannot be inverted 
either, because of a slow l/y 2 -decay of p(y) (see Section 3). 

3. THE DENSITY PROFILE 

A feature of the two-dimensional Coulomb gas is that all coefficients in the low fugacity 
expansion of the pressure and correlation functions are convergent for T > 4 [7]. This 
signals the transition from a conductive phase for T < 4 to an insulator phase for T > 4, in 
the zero-density limit. Similarly, by examining the second moment of the cluster integral 
for the density profile of a single particle, we expect that all the coefficients in the low 
fugacity expansions for the model of subsection 2.1 are convergent for T + 2a > A and 
that this signals the transition from a conductive phase for r + 2a < 4 to an insulator 
phase for r + 2a > 4, in the zero density limit. In this section the low fugacity expansion 
for the density profile will be studied for r + 2a — > 4 + , which is the limit of approaching 
the phase boundary from the dipole side, using the techniques introduced in [2]. More 
precisely, we will study the asymptotic density profile PAcf>(y), which is defined as the 
portion of the low fugacity expansion of p(y) that gives the correct leading order singular 
behaviour of each term in the low fugacity expansion of A0 (2.17) in this limit. 

Alastuey and Cornu [2] complemented their study of the low fugacity expansions of 
the correlations in the dipole phase of the two-dimensional Coulomb gas by an analysis 
of the BGY equations. In subsection 3.5 the first BGY equation is used to compute the 
leading asymptotic behaviour of the density profile. Unlike the resummation of the low 
fugacity expansion calculation, there is no underlying assumption that the phase of the 
model is near the critical point on the dipole side. 

3.1 The expansion at 0(( 2 ) 

Suppose the model of subsection 2.1 is generalized so that each particle is associated with 
a position dependent fugacity ( i— > ((y) (or equivalently is subject to an extra one body 




(2.21) 




(2.22) 
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potential acting perpendicular to the interface). Denote the corresponding N- particle 
canonical partition function by Zn and grand partition function by S. Then from the 
formula 



■*C(v) 



C(»)=C 



it is easy to show that 



Piy) = c m Zl 



+ c 2 



C(y)=C 



6((y) 



«y)=c 1 



+ 0(C 3 ) (3.2) 



C(»)=C; 



Indeed (3.2) applies to any one-component system. Inserting the form of the partition 
functions for the model under consideration we have 

C C 2 f°° f°° i 

Piy)= (2ypiy* + (2y)^y«J-J Xl L ^tMV 
' (x - Xl ) 2 + (y - yi f\ T ' 



x 



-i }+o(c i ) 



(x - xi) 2 + (y + yi) 2 j 
Substituting the term of (3.3) proportional to £, p^\y) say, in (2.17) we find 

C2^ r / 2 rf- r / 2 ~ a + 2 



(3.3) 



A0 (1) /27rg 



T/2 - a + 2 



(3.4) 



(again, and below, we have used the superscript to indicate that only the term proportional 
to this power of ( is being considered). Thus A^ 1 ) is singular in the limit r + 2a — > 4 + , 
and furthermore to leading order is independent of d. Both features are true of Acp^ 
in general. The latter feature implies that only the large-;?/ asymptotic portion of p^ n \y) 
contributes to the leading order singular behaviour of A(f^ n \ and thus p^(y) consists of 
terms in the asymptotic expansion of p^ n \y). With n = 1, there is only one term, which 
is p (1) {y) itself, so trivially p%{y) = p (1 \y). 

Let us now determine the portion of the term of order ( 2 in the low density expansion 
of p(y) which contributes to the leading order singular behavior of A<f>, and thus calculate 

(2) 

PA<t>(y)- For this purpose we consider the double integral which is part of the coefficient of 
the ( 2 term in (3.3) and analyse it for large-y. We first break the integration over y 1 in the 
double integral in (3.3) into two intervals: [d, y] and [y, oo]. A change of variables y\ \— > 
yv 1 shows that the latter interval of integration gives a contribution to the asymptotic 
expansion of p^(y) which is 0(y~ r / 2 ~ a ~( r / 2+a ~ 2 ^), and a corresponding contribution to 
A<// 2 ) which is 0(1/ (r + 2a — 4)). For the interval [d, y], use of the large-?/ expansion 



(x - Xj) 2 + (y - 

(x - xx) 2 + (y + i/i) 2 



r/2 



- 1 



2Tyiy 
x 2 + y 2 



(3.5) 



and integration over x gives a contribution to the asymptotic expansion of p^ 2 \y) of 



c 2 

(2y)ryy 



ry 

/ dy 1 S y (y 1 ) 

J d 



(3.6) 
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where 



s v (yi) ■-- 



2nT 



(2y 1 )r/2yf- 1 ' 
Computing the integral, (3.6) reads 

2vrrC 2 



d<yi<y 



2V /2+a (r/2 + a - 2 



(T/2+0-2) _ y -(r/2+. 



(3.7) 



(3.8) 



Note that the first term in the last brackets above is all that need be included for the 
leading order asymptotic expansion of p^ 2 \y). However both terms give a contribution 
to A<f>W which is 0(1/ (T + 2a- 4) 2 ). Hence p%(y) is given by both terms in (3.8), or 
equivalently the integral formula (3.6). Following [2] we can interpret the integral (3.6) 
as resulting from the partial screening of the fixed particle-image pair of separation 2y by 
the smaller pair of separation 2y 1: via the operator (3.7). 

3.2 Nested dipole hypothesis 

Rather than attempt to calculate p^,(y), n > 3, from the low fugacity expansion (3.2), 
we make a nested dipole chain hypothesis, analogous in idea to that of Kosterlitz and 
Thouless [1] and technically to that of Alastuey and Cornu [2]. Technically we suppose all 
configurations contributing to p^(y) are nested chains of particle- image pairs, with the 
fixed particle- image pair the largest, and the screening operator acting between connected 
particles in the chain only. To specify the chains we can ignore the images and consider 
the different ways of arranging the mobile particles into chains below the root particle 
at y. For example, at 0(( 4 ) there are four distinct configurations, which are illustrated 
graphically in fig. 1. The ordering y > y\ > y2 > y$ > d is equivalent to the nesting of 
the particle-image pairs so that each pair screens a pair of smaller size. The contribution 
to PA^(y) from each graph is 

rv ryi ] r ry ] r ry 

I dy 1 S y (y 1 ) / dy 2 S yi (y 2 ) / dy 1 S y (y 1 ) , / dy 1 S y (y 1 ) 

d J d \ \_J d J \_J d 



yi 



yi 



V2 



, dy 1 S y (yi)[l dy 2 Sy 1 (y 2 )) / dy 1 S y (y 1 ) dy 2 S yi (y 2 ) dy 3 S y2 (y 3 ) 

'd \Jd / J d Jd Jd 

respectively. Furthermore, the graphs need to be weighted by factors of 6,1,3,6 respectively 
to account for relabelling degeneracy, and this linear combination then multiplied by an 
overall factor of C 4 /3!(2y) r / 2 |/ a . 



y - 
yi - 

V2 - 
V3 - 




A 



Figure 1. Graphical representation of the four distinct chains at 0(( 4 
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The nested structure allows the general terms of order n in ( to be calculated by 
recurrence. As is shown in detail in [2, eqs.(4.26)-(4.28)] we have 

&y> = (2y ; /V( „ C _ 1)! ^ 1) W (3*) 



where 



qi+... + qp — n-1-p 

I q (y)= [ V dy'S y (y')S^(y') (3.9c) 



with 



Furthermore, it is shown in [2] that ^ in) as given by this equation can be summed 
over n whatever the form of S y (y'), with the result 



f ry 

pA<p(y) = (2 y y/i y « exp [j d d y^ s y(y^( 2 y^ r yip^ivi) 

Inserting the explicit form (3.7) gives the non-linear integral equation 



(3.10) 



p^(y) = ( 2y )r/2 y a ex P [~ 2nT J d d yiyip^(y 



(3.11) 



which uniquely determines PA4>(y)- Note that even though this equation was derived with 
the assumption r + 2a — > 4+ at each order in (, it is well defined for all values of T + 2a. 

3.3 Evaluation of pA<t>{y) 

By multiplying both sides of (3.11) by y T / 2+a and differentiating with respect to y we 
obtain the non-linear differential equation 

i y 9iy) = -y^=-My)f (3-12a) 

where 

g(y):=y r/2+a p A M (3.12&) 

This is to be solved subject to the initial condition pA</)(d) = (/[(2d) T ^ 2 d a ], obtained by 
substituting y — d in (3. 11). Since the differential equation is first order separable, the 
solution of the initial value problem is straightforward, and we find 

, x c/(2y) r /y 

1 "I" 2 r/2(r/2+a-2)^ a V ' ) 



provided r + 2a ^ 4, while 



c/(2y) r /y 

= l + ^-r/^rClogy (3 - 136) 



for T + 2a = 4. 
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Inspection of (3.13a) show that the exponents in the y-dependent terms of p^(y) 
depend on T/2 + a only. In particular there is no dependence on the fugacity ( and 
consequently the phase transition will be independent of (. This behaviour is to contrasted 
with the two-dimensional Coulomb gas, in which the powers in the decay of the asymptotic 
charge-charge correlation at order (the quantity analagous to pA<f>(y)) are dependent 
on ( and the phase transition is ^-dependent. 

We can use the resummations (3.13) to calculate A(j) to all orders in (. First, from 
(2.17) we see that we require 

PA<t>{y) = e > 

for A0 to be finite. Now from (3.13) we have 

2 ~l^ a , T + 2a<4 
PaAv) ~ { r + 2« = 4 (3.14) 



y r/2+a ■ 



T + 2a > 4 



Hence A0 is finite for T + 2a > 4 (the dipole phase), and infinite for T + 2a < 4 (the 
conductor phase), independent of the fugacity ( as anticipated above. 

The explicit value of A(p in the insulator phase is given by computing the integral 
(2.17) with p(y) replaced by (3.13a). We find 



A0/27rg= ^ log 



2nT(d 2 ~ r / 2 - a 
1 + 2 r / 2 (r/2 + a - 2) 



(3.15) 



Again we emphasize that even though the intermediate steps leading to this result require 
r + 2a — >• 4+ at each order in (, (3.15) is well defined for all r + 2a > 4. 

3.4 Renormalization flow equation 

In the two-dimensional Coulomb gas the renormalization flow equation relates the asymp- 
totic charge-charge correlation to the space dependent dielectric constant, with the space 
variable an implicit parameter. In the present model we can obtain a renormalization 
flow equation by relating pA<t>(y) to 

ry 

A(j)(y) := 2irq / dy 1 y 1 p A 4>(yi) (3.16) 

J d 



To do this we differentiate (3.16) and use the definition (3.12b) to obtain 

2nqy 1 - r / 2 - a g(y) (3.17) 



d^iv) _ o. L- 1.72-n 



dy 

Dividing (3.12a) and (3.17) then gives the desired equation 



My) - T - a (v) (3.18) 



dA(j)(y) q- 

This flow equation is subject to the initial condition g(d) = 2~ r / 2 ( and A<j>(d) = 0, and 
its exact solution is thus 

g( y ) = 2 -r/2 Ce -rA^)/ 9 ( 319 ) 
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From (3.15), for T+2a > 4 the allowed values of A<f>(y) is the finite interval [0, A(f>/2nq}, 
while for T + 2a < 4 A(f)(y) takes on all values in [0, oo[. The flow diagram obtained from 
(3.19) thus has the appearance sketched in figure 2. 



Figure 2. The flow diagram, where we have written x(y) := TA(y)/q. The different trajectories 
correspond to different values of T + 2a. The trajectories terminate for T + 2a > 4. 

3.5 The first BGY equation 

In this subsection we will complement the above low fugacity resummation study by an 
asymptotic analysis of the first BGY equation. Let us denote the force on a particle at f[ 
due to a particle at f* 2 by F 2 ±, so that 



'21 



-Vi0(ri,r£) (3.20) 



where 4>(f{, 7%) is given by (2.1a). Furthermore denote the force on a particle at f[ due to 
the self image particle and the one body potential by F™ and Ff xt respectively so that 

2 

F im = and ^cxt = ■ (321) 

Then in terms of these forces the first BGY equation for the system is 

Vip(rl) = PFf*p{f[) + f3F~i m p(fl) + (3 [ dr 2 F 21 p {2) (rl, r 2 ) (3.22a) 

Let us consider the y-components of this equation for y\ — > oo. We might expect that 
in this limit we can replace p^(fi,f 2 ) in the final term in (3.22a) by p(fi)p(f 2 ), which is 
equivalent to saying that if we write the final term in (3.22a) as 



Jv 



dr 2 F 21 p{n)p{f 2 ) + (3 [ df 2 F 2 ip (2)T (n, r 2 ) (3.226) 

Jv 



the term involving p^ 7 '(fi, r 2 ) decays faster for large- y± than the term involving p(fi)p(f 2 ) (which 
is a mean-field term), plus the one body forces on the r.h.s. of (3.22a) . In the appendix 
this latter statement is proved subject to a mild clustering assumption. Thus, neglecting 
the second term in (3.22b), we are left with the mean-field type equation 



T + 2a 

+ 



2yi 



P [ dr 2 (F 21 ) y p(y 2 ) 
Jv 



p(Vi) (3-23) 
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The integral over T> in this equation can be simplified: 

/ dr 2 (F 2l ) y p(y 2 ) 
Jv 



-Q 



-q 



-oo 
oo 



dy 2 



^vMn-r 2 \) + ^v c (\n-r2\) 



P{V2) 



dy 2 

= -2vrg 2 / dy 2 p(y 



d d 
— v c (0; yi - y 2 ) - t—v c (0; y x + y 2 ) 
dyi dyi 



(this result can also be derived from Gauss's theorem in electrostatics). The mean-field 
equation now reads 



d_ 

dyi 



P(yi) = 



T + 2a 



2irT 



dy 2 p(y2 



p(yi 



To solve this equation for large-yx we seek a solution of the form 



Substituting this in (3.24) gives 

pc 



p(Vi) - — p 



[V + 2a)c 



(3.24) 



(3.25) 



27rrc 2 



P+i 
Vi 



2y p i +1 (p - i)vi 



2p-l 



(3.26) 



For p > 2 the second term on the right hand side of (3.26) can be ignored and we obtain 
a solution provided 



p= {T + 2a)/2 (and thus T + 2a > 4) 

For p — 2, (3.25) is an exact solution of (3.24) provided 

4 - T - 2a 



c = 



and 



T + 2a < 4 



(3.27a) 



(3.276) 



For r + 2a = 4 we find, after equating the first two orders on both sides of (3.24), that 

1 



p(yi) 



2ixTy 1 2 \ogy l 



(3.27c) 



is an asymptotic solution. 

We emphasize that the above analysis is asymptotically exact and non-perturbative: 
the asymptotic formulas obtained hold for all values of T, a and (. As such, we can use 
these results to test the predictions (3.14) for the leading asymptotics of p(y) as derived 
from pA<t>{y)- Surprisingly the results obtained from p^{y) are in complete agreement 
with the non-perturbative exact results, even though it has been assumed in the derivation 
of PAcf>(y) that the phase of the model is near the zero-fugacity critical point on the dipole 
side. 
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4. Comparison with the solvable case 



4.1 The phase 

When T = 2 the model of subsection 2.2 is exactly solvable for all a [5]. The exact 
expressions for the density profile and truncated two particle distribution function are 

roc e -4iryt 



and 

e 2nixt g-27r(yi+j/ 2 )t 



p 1 (f 1 ,f 2 ) = -(2 7 rC) 2 (y 1 y 2 ) 



eft- 



(4.2) 



i + 2ir( dyy-«e- w ' 

These expressions were used in [5] to show that for a < 1 the dipole moment of the 
internal screening cloud D(y ) as defined by (2.18) vanishes, while for a > 1 it is non- 
zero. This behaviour was interpreted as indicating that the system exhibits a conductive 
phase for a < 1 and an insulator phase for a > 1. In subsection 2.2 we have showed that 
the true indicator of a conductive phase is the sum rule (2.15), and the vanishing of D(y ) 
in a conductive phase is a corollary of this stronger requirement. 
Noting that for a one-component system 

S(y, y'- x-x') = q 2 [p(y f )5(x - x')S(y - y') + p T (y, y'- x - x')] (4.3) 
and using the exact results (4.1) and (4.2), it is a straightforward exercise to show 

poo roo roo 1 _|_ Oirf i' V0 HYY~ a 

*»/. *L"S M ;*- = 1- + ' „ (4.4) 



For ct < 1 the second term on the r.h.s. of (4.4) vanishes and the sum rule (2.15) holds, 
thus implying a conductive phase. For a > 1 (2.15) is not obeyed so the phase in an 
insulator. These conclusions are in agreement with those reached in [5]. 

In subsection 3.3 we have shown that the potential drop A0 diverges for r + 2a < 4 
but is finite for r + 2a; > 4. In accordance with the interpretation of the formula (2.17) as 
saying A0 is proportional to the mean distance of separation within the dipole formed by 
a particle and its image, we have taken this behaviour to be an alternative phase indicator 
to the sum rule (2.15) for this system. For the solvable model we have the exact result [5] 

A0/27rg=-^log(l + 27rC^4), a > 1 (4.5) 
47r y a — 1 / 

Remarkably, the expression for A<p (3.15) with T = 2 deduced from the asymptotic density 
profile PA<f>(y) is in precise agreement with this exact expression. 

4.2 The asymptotic density pA<p(y) 

The asymptotic density pA<p{y) is defined as the portion of the asymptotic expansion of 
p(y) that gives the correct singular behaviour of A<p as T + 2a — > 4 + at each order in (. 
From (4.1) we can calculate pA<f>(y) exactly at T = 2. 

Expanding (4.1) as a power series in ( and then performing the integration over t gives 

p(y) - ^g(- w*+'> S WyM , % . ■ ■ f + . . . + v, + (4.6) 
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For large-y the final factor in the integral can be approximated by 1 and we obtain 



Zy 3=0 



/■ oo 

V^/ r,_A\1..-j(a+l) 

\ cc — 1 y 

(4-7) 

The correction term in the asymptotic expansion above does not contribute to PA<j>(y)- 
Ignoring this term, we see that a geometric series remains, which after summation gives 

Comparison of this exact result at T = 2 with the result (3.13a) obtained from the low fu- 
gacity resummation using the nested dipole chain hypothesis shows that the resummation 
is exact at this coupling. 

4.3 Leading asymptotics of density profile 

As noted in [5], the leading large-y behaviour of the density profile at T = 2 is readily 
computed from (4.1). We find precise agreement with the behaviour (3.14) at T = 2, 
which is obtained from both the low fugacity resummation and the mean-field equation. 

It is interesting to note that since the asymptotic form of the density profile in the 
conductive phase is 

. , 4-T-2a . . 

p(y) ~ -Siy- (49) 

the phase transition occurs when the 1/y 2 tail vanishes. 



5. CONCLUSION 

The metal wall one-component plasma model of subsection 2.1 exhibits both a conductive 
and insulating phase. It has the special property of admitting an exact solution for the 
thermodynamics and all correlations along the line (T, a) = (2, a) in parameter space [5]. 
This line intersects the transition line r + 2a = 4. For general values of the parameters 
the transition can be analysed in a similar way to the Kosterlitz-Thouless transition in 
the two-dimensional Coulomb gas [2]. In particular, by making an hypothesis that the 
dominant configurations are nested dipole chains (this is checked explicitly at 0(( 2 )), the 
low fugacity expansion of the asymptotic density PA<p(y) can be resummed, and calculated 
explicitly as the solution of a non-linear differential equation. 

Comparison with exact solution verifies that the general expression for pA<f>(y) is ex- 
act at T = 2. This provides compelling evidence for the correctness of the underlying 
nested dipole chain hypothesis. Since the nested dipole chain hypothesis also underlies 
the iterated mean-field equations of Kosterlitz and Thouless [1] (which are equivalent to 
the Kosterlitz renormalization equations [8]), we have also added further weight to the 
validity of these equations. 
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Appendix 

In this appendix, we will prove that in the limit y\ — > oo the second term in (3.22b), 

P{ 2) := (3 dr 2 F 21 pJ 2) (n, r 2 ) (Al) 

decays faster than the sum of the first term and the one body terms on the r.h.s. of the 
BGY equation (3.22a), and therefore can be neglected in this limit. Our analysis is based 
on the simple assumption that for some 1 > e > 

|pf 2) (rl,f 2 )| < p(rl)p(r 2 ) (r^r) (A2) 

\\n - r 2 \J 

where / is a given length. We stress that the hypothesis (A2) is very reasonable since it 
merely asserts that the Ursell function pj 2 ^(fi, f 2 ) / [p(fip(f 2 )} decays for large separations 
\f\ — r 2 \ at least as fast as an inverse power. This weak clustering property surely holds 
in any homogeneous or inhomogeneous fluid phase. 

Using (A2) and the inequality \y 2 — y\\ < \r 2 — fi|, we find 



|^ (2) | 



<2r P (y 1 )J v df 2P (y 2 )(^-^ ) ) y 2[( , 2 _ :Ci)2 + (y2 _ yi)2]1/2[(:C2 _ :Cl)2 + (y2 + yi)2] V 2 

(A3) 

In the integral on the r.h.s. of (A3), we can perform the integration over x 2 according to 



/: 



^oo X2 i( X2 - Xl )2 + (y 2 - y^Kx, - Xl )2 + (y 2 + (y 2 + y ± ) [y 2 + y ± ) 

(A4) 

where K(k) is the complete elliptic integral of the first kind 

*/2 dd 



A/1-A; 2 sin^ 



K(k) 

By splitting the domain of integration over y 2 into the intervals [d,yi/2] and [yi/2,oo[, 
we then find from (A3) 



+4Tp( yi ) r dy 2 p(y 2 ) ( *—) * K (h/Ml) (A5 ) 

Jyi/z \\y2-y1\J \y2 + yij 

(we have also used the monotonicity of K{k)). 

For yi large, p(yi) is expected to decay as c/y[, with p > 2 (this is shown explicitly in 
subection 3.5). The integral 

/ dy 2 y 2 p(y 2 ) 

J d 
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then remains bounded by some constant times log(yi). Therefore the first term on the 
r.h.s. of (A5) decays at least as fast as \og(yi)/yi +1+e which is faster than the one-body 
self image and external potential terms (~ l/y r { +1 ) appearing on the r.h.s. of the BGY 
equation (3.22a). Also, the integral 



remains bounded by a constant times 1/yf t , as shown by the variable change yi = ay\. 
Indeed, the dimensionless integral 



is finite because the singularity of K(2a 1 ^ 2 / (a + 1)) at a — 1 is only logarithmic: 

K(2a 1/2 /(a + 1)) log|a- 1| 

when a — > 1. Then, the second term on the r.h.s. of (A5) decays at least as fast as 
l/y\ p ~ l+e , which is faster than the decay of the first term in (3.22b) (~ l/yl^ 1 ). Thus 
the whole two-body force (Al) can be neglected with respect to the other terms of the 
BGY equation (3.22a) in the limit y\ — > oo. 





17 



REFERENCES 



[1] J.M. Kosterlitz and D.J. Thouless, J. Phys. C 6 (1973) 1181 
[2] A. Alastuey and F. Cornu, J. Stat. Phys. 66 (1992) 165 

[3] J. Clerouin and J.P Hansen, Phys. Rev. Lett. 54 (1985) 2277; A. Alastuey, F. Cornu 
and B. Jancovici, Phys. Rev. A 38 (1988) 4916 

[4] A. Alastuey, Mol. Phys. 52 (1984) 637 

[5] P.J. Forrester, Int. J. Mod. Phys. A 7, suppl. 1A, (1992) 303 

[6] B. Jancovici, J. Phys. (Paris) 47 (1986) 389; Ph. A. Martin, Rev. Mod. Phys. 60 
(1988) 1075 

[7] E.R. Speer, J. Stat. Phys. 42 (1986) 895 
[8] A.P Young, J. Phys. C 11 (1978) L453 

[9] I.M. Ryshik and I.S. Gradstein, "Tables of Series, Products and Integrals" (VEB 
Deutscher Verlag der Wissenschaften, Berlin, 1963) 



18 



